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Phase coexistence in polydisperse mixture of hard-sphere colloidal 
and flexible chain particles 
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Abstract. - A theoretical scheme for the calculation of the full phase diagram (including cloud 
and shadow curves, binodals and distribution functions of the coexisting phases) for colloid- 
polymer mixtures with polymer chain length polydispersity and hard-sphere colloidal and poly- 
meric monomer sizes polydispersity is proposed. The scheme combines thermodynamic pertur- 
bation theory for associating fluids and recently developed method used to determine the phase 
diagram of polydisperse spherical shape colloidal fluids (L.Bellier-Castella et al, J.C'hem.Phys. 
113, 8337(2000)). By way of illustration we present and discuss the full phase diagram for the 
mixture with polydispersity in the size of the hard-sphere colloidal particles. 



Introduction. — Since the pioneering studies of 
Asakura and Oosawa [1,2] and Vrij [3] substantial amount 
of efforts have been focused on the development of the 
theoretical methods describing the phase behavior of the 
athermal colloid-polymer mixtures (see Refs. [4-6] and 
references therein). Usually in the vast majority of the 
theories developed so far interaction between polymers 
is either ignored or treated using approximations differ- 
ent from those assumed to describe polymer-colloid and 
colloid-colloid interactions. This feature imposes certain 
restrictions on the possibilities of the theory, for exam- 
ple the theories, which ignore polymer-polymer interac- 
tions are restricted to the mixtures of short polymers 
and large colloidal particles (ideal polymer limit). Re- 
cently an attempt to describe the phase behavior of the 
colloid-polymer mixture with both components treated on 
an equal footing using polymer reference interaction site 
model integral-equation approach [7] and thermodynamic 
perturbation theory (TPT) for associating fluids [6] was 
made. In these studies polymer molecules were modeled 
as a flexible chains of tangentially bonded hard-sphere 
monomers and colloidal particles were represented as a 
hard spheres. Later TPT approach was extended [8] to 
account for the chain length polydispersity effects. In 
this study we propose further extension of the TPT to 
account for polydispersity in the polymer chain length 



and in the hard-sphere sizes of both colloidal particles and 
polymer monomers. This extension is based on the previ- 
ously developed method [9, 10], which allows us to deter- 
mine the full phase diagram (including binodals and cloud 
and shadow curves) and to discuss fractionation effects on 
the level of the distribution functions of the two daughter 
phases. 



The model. — We consider polydisperse mixture of 
hard-sphere flexible chain particles represented by m tan- 
gentially bonded hard spheres of diameter a. In the fol- 
lowing we will distinguish between p— type of the particles 
(polymers) and c— type of the particles (colloids). Thus 
the species of each particle is characterized by the set of 
three variables (a, m, a) with a denoting the type of the 
particle (either p or c), m = 1,2, oo and < a < oo. 
The number density of the a— type of the particles is p a 
and the overall number density p is p = p p + p c . The 
species variables a, m and a are distributed according to 
the distribution function F a (m,a), which is positive and 
satisfies the following normalizing conditions 



d ° F a.{m,<T) = I. 
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Further we put F a (m,a) = a a f a (m,a), where a a denote 
the fraction of the a— type of the particles 
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obviously, a p + a c = 1 and partial distribution functions 
f a (m,a) are normalized. 

Thermodynamical properties. — Thcrmodynami- 
cal properties of the model at hand are calculated using 
TPT of Werthcim [11, 12]. According to TPT Helmholtz 
free energy of the system A is represented as a sum of 
three terms A = A ic i + Ah s + A c h, where An is the ideal 
gas term 



Chemical potential is obtained as the functional derivative 
of the free energy density F/V with respect to the function 
p a {m,a) = p a F a (m, a): 

fi a {m, a) = In [p a F(m, a)} + p a ,hs( m , cr ) + Va,ch(m, cr), 
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and A c /j is the term describing formation of the chains 

a m 
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Here V is the system volume, k is the Boltzmann con- 
stant, T is the temperature, Co>Ci)C2 5 C3 are the distribu- 
tion function moments 

C™ = ^PJ2J2 m J da F a(m,cr)a n , n = 0,l,2,3; 



A = 1 — (3, gda S \&) is the hard-sphere contact value 
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In (4) and (7) Boublik's [13] and Mansoori's et al. [14] ex- 
pressions for Ahs an d gia^ {&) have been utilized. All the 
rest of thermodynamical properties can be obtained from 
Helmholtz free energy using the standard thermodynam- 
ical relations. Taking the volume derivative of the free 
energy we get the following expression for the pressure: 
PP = P + PP hs + pPch , where 
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One can easily see that thermodynamical properties of the 
model at hand are defined by the set of a finite number of 
the distribution function moments, i.e. four regular mo- 
ments (0X1X2,(3 an d two generalized moments (O, *). 
Thus polydisperse mixture of the chain particles treated 
within TPT belong to the class of truncatable free energy 
models [15]. 

Phase equilibrium conditions. We assume that 
at a certain density p^ and composition F^ (to, a) the 
system separates into two phases with the densities p^> 
and p( 2 \ and compositions F^[m, a) and F^ (to, a). 
Hereafter the upper index (0) refer to the parent phase 
and the upper indices (1) and (2) refer to the daughter 
phases. At equilibrium these quantities take the values, 
which follows from the phase equilibrium conditions, i.e.: 
(i) conservation of the total volume of the system, (ii) 
conservation of the total number of the particles of each 
species, (hi) equality of the chemical potentials of particles 
of the same species in the coexisting phases, (iv) equality 
of the pressure in the coexisting phases. These conditions 
finally lead to the following set of relations [9, 10]: 



F^(m,a)=FW ( m ,a) 
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^2^2 j da F^ a \m,a) = 1, for a = 1 or a = 2, (19) 
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a = 1, 2. We have 




Fig. 1: Phase diagram in: (a) P* vs p* and (b) P* vs a p 
coordinate frames. Cloud and shadow curves are represented 
by thick and thin solid lines, respectively, binodal curves with: 
(1) o4 0) = 0.541 and (2) a p 0) = 0.931 by the broken lines and 
with a p = a p °l r = 0.744 by the dotted line. The dashed- 
dotted line denotes the binodal curve for the bidisperse version 
of the model and critical points are indicated by the diamonds. 
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/ii"'' is the excess (over the ideal gas) chemical potential 
of the particle (a, m, a) in the phase a and [. . .] denote 
functional dependence. The relation between (m, a) 
and daughter phase distribution function Fa a \m, a), i.e., 
Eq. (17), follows from the phase equilibrium conditions 
(i)-(iii). 

Relations (17)-(19) represent a closed set of equations to 
be solved for the unknowns p(°^ and Fa (to, a); this set 
have to be solved for every value of the species variables 
a, to and a. However, since thermodynamical properties 
of the model at hand arc defined by the finite number 
of the moments we can map this set of equations onto 
a closed set of fourteen algebraic equations for p^ and 
twelve moments (n — 0, 1, 2, 3), f^ Q ) and ^> (a \ where 
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a m 
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xQW^mW;^ 1 )},}^)}), (24) 
where {X( Q )} represent unknowns of the problem, i.e. 

{lW} = { / )( a ) 1 C,fl (a) 1 * (a) }, n = 0,1,2,3; a = 1,2. 

The remaining two equations follows from the equality of 
the pressure in coexisting phases (18), 

(>); [X^Y) = p(2) (V 2) ; {^ (2) }) . (25) 



and from the normalizing condition (19) for either phase 
a = 1 or a = 2, 



a m 
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(26) 



Solution of the set of equations (22)-(26) for a given den- 
sity and distribution function F^ (to, a) of the par- 
ent phase gives the densities and distribution func- 
tions Fa a \m, a) of the two coexisting daughter phases. 
The coexisting densities at different densities of the parent 
phase p(°) defined binodals, which are terminated when 
a density of one of the phases is equal to the parent 
phase density p^ . For the different values of the colloidal 
(and polymer) particles fraction (a^ = 1 — a^) in 
the parent phase these termination points form the cloud 
and shadow coexisting curves. These curves intersect at 
the critical point, which is characterized by the critical 



density p c 



,(i) 



,(2) 



//°) and critical colloidal 



(and polymeric) composition a C;Cr = a^P = = 
(ct p ,cr = 1 — Oi CiCr ). The cloud-shadow curves can be ob- 
tained as a special solution of the general coexisting prob- 
lem, when the properties of one phase are equal to the 
properties of the parent phase: assuming that the phase 
a = 2 is the cloud phase, i.e. p^ = p(°\ and following the 
above scheme we will end up with the same set of equa- 
tions (22)-(26), but with p^ and Pi 2) (TO,cr) substituted 
by and Fa°\m,a), respectively. 
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Fig. 2: The same as in Figure 1 in: (a) r\ vs p* and (b) r\ v 
vs r\c coordinate frames. Thin solid tie lines and broken tie 
lines connect the points on the cloud and shadow curves and 
dotted tie lines connect the points on the binodal curve for the 
bidisperse version of the model. 



Results and discussion. — To illustrate the scheme 
developed in the previous sections we present here numer- 
ical results for the phase behavior of the polydisperse mix- 
ture of colloidal and polymeric particles with poydispersity 
in the size of the colloidal particles only. Colloids are rep- 
resented by the polydisperse mixture of hard spheres and 
polymers are modeled by the hard-sphere flexible chains 
with fixed chain length mo and hard-sphere size of the 
monomeric units a p . Thus for the distribution function of 



the parent phase Fa (m, a) we have 

F^(m, a) = a ( p ^S (a - cr p ) 8 m ^ ma 



(27) 
(28) 



where for the colloidal diameter distribution fc°\cr) we 
have chosen beta-distribution, given by 



f?\a)=B-\ % v) U 
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where 0(x) is the Heaviside step function, B(-f, v) is the 
Beta function [10] and 7 and v are related to the first 



Fig. 3: Distribution functions of the parent phase (dotted line), 
low density phase (solid line) and high density phase (dashed 
line) on the shadow curve at P* = 8. Here a* — a/(a}^ 
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(29) 

Calculations were carried out for the model parameters 
chosen to be: m = 50, a p /(a)^ = 0.06, D a = 0.1 and 
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At sufficiently low density of the parent phase 
the mixture is stable as a single phase. As is in- 
creased the system will phase separate into two phases: 
low density phase and high density phase. In Figures 1 
and 2 we show the phase diagram in a different coordi- 
nate frames, i.e. P* vs p* (Fig. la), P* vs a p (Fig. 
lb), 77 vs p* (Fig. 2a) and r\ p vs i] c (Fig. 2b) , where 

P* = /3P((cr)c ^) 3 , p* — p((o-}c°') 3 , V and r) a are the over- 
all packing fraction and packing fraction of the a— type of 
the particles, respectively. The phase diagram shown in 
Figure 1 includes the cloud and shadow curves, along with 
three binodals for three selected parent phase fractions 



„(°) 



= 0.541, 0.744, 0.931, one of them being the criti- 



cal fraction a 



(o) 

p,cr 



0.744. The two branches of the bin- 
Jo) 



odals for the two selected a p '' values (cip ^ = 0.541, 0.931) 
terminate on the cloud and shadow curve, where the p 
values of the end points of the respective binodals on 
the cloud curve are equal to the p^ values of the bin- 
odals in these points. For the critical value of the par- 
ent phase fraction a p °l r = 0.744 corresponding binodals 
meet in the critical point. In Figure 2 we include only 
the cloud and shadow curves. For the reference in both 
figures we also display the binodal curve for the bidisperse 
version of the model (D^ = 0). Comparison between 
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Fig. 4: Average hard-sphere size a ( c a) = {a) ( ^/{a}^ (a) and 
distribution function width L>S (b) of the colloidal particles 
along the shadow curve. Solid lines represent low density phase 
and dashed lines denote high density phase. Critical points are 
indicated by the diamonds. 



the bidisperse and polydisperse versions of the colloid- 
polymer mixture shows that polydispersity extends the 
region of the phase instability, shifting the critical point 
to the lower values of the pressure and density (P* r bid 



ments, (cr)^ and (29). These two quantities char- 

acterize size-distribution of the particles and distribution 
function width, respectively. In Figure 4 we present (cr)^ 
and along the shadow curve. It can be seen that the 
larger size particles prefer the lower density phase and 
the smaller size particles arc predominantly encountered 
in the higher density phase. As the pressure increases we 
observe increase in the mean size of the particles and de- 
crease in the distribution function width in the low-density 
phase and decrease in the mean size of the particles and 
increase in the distribution function width in the high- 
density phase. The magnitude of these changes for the 
mean size of the particles in both phases are similar (Fig 
4a): we observe a strong sharpening of the low-density 
phase distribution function with increasing pressure, while 
increase in the width of the high-density phase distribu- 
tion function do not exceed 13% of its parent phase value 
(Fig. 4b). A more detailed and systematic investigation 
of these effects for the mixture with polydispersity in both 
colloidal and polymeric subsystems will be postponed to 
a future contributions. 

Conclusions. — With the concept presented above 
we are able to calculate the full phase diagram of poly- 
disperse hard-sphere colloidal and flexible chain particles 
with all types of polydispersity. For a particular system 
with polydispersity in the size of the colloidal particles we 
present and discuss the phase diagram and corresponding 
fractionation effects. 
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(bid) 



0.776, a 



(polyd) 
p,cr 
(bid) 



0.744, Fig. lb) and poly- 
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tion (a 

mcr packing fraction (rj p "c^' = 0.009, 
2b). Corresponding relative changes of the colloidal (rj CtCr ) 
and overall (r] cr ) critical packing fraction are negligible. 

According to Figure 2 packing fraction of the low den- 
sity phase is higher than that of the high density phase; 
thus the density of the liquid-like phase is lower than that 
of the gas-like phase. This is not surprising as the frac- 
tion of the colloidal particles in the low density phase is 
larger than in the high density phase. In Figure 3 we 
show daughter phase distribution functions fc a \<?) on the 
shadow curve at P* — 8 together with the parent phase 
distribution function /i° (<r). More detailed information 
about the composition of the coexisting phases can be ex- 
tracted from the analysis of the distribution functions of 
the two daughter phases in terms of their first two mo- 
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